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Overview 

• Part I: Study Objectives & Approach 

• Part II: Background Info 

• Water Chemistry 

• Reactive Transport Model (OTEQ) 

• Part III: Results 

• Concentrations and WQ Standards 

• Loads and Sources 

• OTEQ calibration and remedial scenarios 
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Part I: Study Objectives & Approach 
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Study Objective 

Main Objective: 

Use a reactive transport model to evaluate 
remedial options for Cement Creek 

Data Requirements: 

- Synoptic data set representing “steady-state” conditions 

- Spatial profiles of streamflow and concentration 

October 2012 Synoptic 



Preliminary Results as of 
September 30, 2014 



1805684 


ED 000792 00006514-00004 




Approach 

Develop spatial profiles of streamflow & concentration 

Study Reach: Ross Basin to A72 (1 1 .5 miles) 

October 2012: i l|||||l||||l|f|| 

- Acoustic Doppler Velocimeter (ADV) to estimate streamflow 

- 2 sampling teams & 3 ADVs 

- leapfrog sampling and day-to-day overlap (replication) 

Loading analysis & modeling assume streamflow and 

concentration do not vary with time. 



Preliminary Results as of 
September 30, 2014 
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October 2012 Sample Locations 


Chemical “snapshot” 

- 30 stream sites 

- 15 inflows 

Sampled for: 

- pH, Alkalinity 

- cations (Fe, Al, Zn...) 

- anions (S0 4 , Cl...) 

leapfrog sampling 
day-to-day replication 

flow from ADV 


Red & Bonita 


American Tunnel 


Preliminary Results as of 
September 30, 2014 
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Porcupine Gulch 


linary Results as of 
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October 2012 
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Part II: Background Info: 

Water Chemistry & Reactive Transport Model (OTEQ) 
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Cement Creek, October 201 


Reactive Transport 


Reaction / Removal 
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iliminary Results as of 
ptember 30, 2014 


.5 

5.5 

6.5 ?? 


pH 



1805684 


ED 000792 00006514-00012 













pH, October 2012 


Reactive Transport 
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Preliminary Results as of 
September 30, 2014 
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Preliminary Results as of 
September 30, 2014 
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Reactive T ransport 
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One-dimensional Transport w/ EQuilibrium Chemistry 


Transport (OTIS) 


Advection + Dispersion + Inflow + Storage 


Equilibrium Chemistry (MINTEQ) 


http://water.usgs.gov/software/OTEQ 
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linary Results as of 
ber30, 2014 
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Results 
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October 2012: Data Set Quality Assurance 

• Loading analysis & modeling assume streamflow and 
concentration do not vary with time (Steady State) 

• Data “Backbone” - Spatial Streamflow Profile 

Essential starting point for loading analysis & modeling 

- ■ leapfrog approach indicates streamflow was -Steady 

- r. means of estimating streamflow provio ' redundancy; 
mass balance calculations used to develop spatial flow profile. 



• Multiple observations of concentration provided by 
leapfrog approach also indicate steady-state 


Zn Concentration on different days 

: cheek for day-to-day cone, variation 


: □ ADV measurements 
— hi Line 


C2 > ("I (cone, increase over time) j 


( l > C2 tconc. decrease over time) ; 


Preliminary Results as of 
mber 30, 2014 
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Part III: Results 

Concentrations and Water Quality Standards 
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Part III: Results - Loads & Sources 


Synoptic Study: 

Spatial profiles of Streamflow & Cone. 

Load (mass / time) = flow * concentration 

— ► Spatial profiles of mass load 

+ changes in load used to identify 
source areas 


OUTLET 

Synoptic Sampling 
October 2-4, 2012 
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Part IV: Results - Loads & Sources 
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Results - Loads & Sources 


Percent Contribution to S04 load at A72 

October 2013 {19 April 2013 Q) 


Percent Contribution to Zn load at A72 

October 2012 (19Apr2013 Q profile) 


Stream Segments Bracketing Source Areas 


Stream Segments Bracketing Source Areas 
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Part III: Results - Loads & Sources 
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Oct. 2012: Sources of Zinc 
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Oct. 2012: Sources of Iron 
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Oct. 2012: Sources of Cadmium 



Preliminary Results as of 
September 30, 2014 
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Oct. 2012: Sources of Aluminum 
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Part III Results: 

Reactive transport modeling 
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Study Objective 
Main Objective: 

Use a reactive transport model to evaluate 
remedial options for Cement Creek 

Data Requirements: 

- Synoptic data set representing “steady-state” conditions 

- Spatial profiles of streamflow and concentration 

Doing the 2-Step: 

- Model Calibration: Reproduce existing conditions 

- Prediction: Modify calibrated model to reflect remediation 


111 


Preliminary Results as of 
September 30, 2014 
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Model Calibration: Reproduce existing conditions 


Components: 

Al, As, Cd, Cu, Fe(ll), Fe(lll), H+, Ni, Pb, S0 4 , Zn, Ca, CO : 

Reactions: 

precipitation: Fe(0FI)3, AI(OFi) 3 

sorption of As, Cd, Cu, Pb, and Zn onto FIFO 

oxidation of Fe(II) 

degassing of C02 

Model Input: 

spatial streamflow profile (1 flow / site based on QA) 


inflow chemistry (synoptic sampling) 
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Model Calibration: Inflow Chemistry (C £ ) 

Set C L to reproduce stream data; Zinc input values realistic 


Zinc Concentrations — Stream & Inflow 

October 2012 
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Model Calibration: Modify Reaction Constants 

adjust AI(OH)3 Ksp to reproduce precipitation below Cement/Animas confluence 


A1 Calibration 


Total Recoverable (unfiltered) 
Dissolved (unfiltered) 

Tot Recov (no degas calibr) 
Tot Recov (degas calibr) 
Dissolved (no degas calibr) 
Dissolved (degas calibr) 


Adjust K s to 
reproduce precipitation 


'eliminary Re: 
leptember 30 


1805684 


ED 000792 00006514-00033 



Predicting Effects of Remediation 

• Modify calibrated model to reflect proposed 
remedial action 

• Remedial Scenarios: 

- Plug Red & Bonita Adit (100% reduction) 

- 1996 Treatment Plant Conditions 





Predicting Effects of Remediation 

Role of Uncertainty 


Sources of Uncertainty: 

- Equilibrium “constants” describing precip tati n 
and sorption vary over a wide range 

- degassing of C0 2 may or may not be an 
important process effecting pH 

- variation in low-flow water quality 


ill 



Role of Uncertainty 
Factor #1 : Equilibrium Constants 

• Solubility Products for Precipitation of Al & Fe 

- known to vary over a wide range 

- default values modified during calibration 

• Surface complexation constants for sorption of 
As, Cd, Cu, Ni, Pb, Zn 

- defaults based on best estimates of Dzombak & Morel 

- default values modified during calibration 

• Dealing w/ Uncertainty: 

- consider predictive simulations w / both 
calibrated values and defaults 


Role of Uncertainty 
Factor #2: degassing effect on pH 


• pH can be difficult to measure; pH data are uncertain 

• Calibrated models with and without degassing have 
been developed 

• Dealing w / Uncertainty: 

- consider predictive simulations w/ and w/o degassing 



Dissolved Zn [mg/L] 



Role of Uncertainty 

Factor #3: Variation in low-flow water quality 

Inflow concentrations are set during calibration 
Inflow concentrations are known to vary: 

Animas River above Cement (A68) Dealiliq W/ Uncertainty: 


Conduct “worst case” 
simulations where Mineral 
Creek and Upper Animas 
inflow concentrations are set 
equal to the max 
concentrations ( 2006 / 2012 ) 
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Dealing w/ Uncertainty 

• Sensitivity Analysis 

• Equilibrium constants: calibrated vs. defaults 

• Degassing: with vs. without 

• Inflow concentrations: calibrated vs. worst case 

• 8 simulations for each remedial scenario 

- provides range of possible outcomes 

- “best estimate” = calibrated degas model (B) 

- compare results based on median (M) 




Remedial Scenario #1 : 100% removal of R&B load 


Modify the calibrated model 
such that no water enters 
between the two stream sites 
that bracket the Red & Bonita 
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Results: 1 00% removal of Red & Bonita load 


Cd, 100% Removal of Red & Bonita load 


Fe, 100% Removal of Red & Bonita load 
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Results: 1 00% removal of Red & Bonita load 


100% Removal of Red & Bonita Load 

Dissolved Zinc at A72 


Observed data 


Chronic Standard 



Range 

M 

Median Simulation 

B 

Best Estimate Simulation 


1805684 


ED 000792 00006514-00042 



Normalized Results 

• Concentrations vary between constituents 
- normalize to allow for comparison 

C n = 100 * (Prediction - Standard) 

(Observed Data - Standard) 

If: 

Prediction = Observed Data (2012) — > C n = 100 
Prediction = Standard —> C n = 0 




Normalized Results: 

100 % Removal of Red & Bonita load 


Existing conditions 


Standard 


Note: Al and Fe results assume all precipitated mass is removed (standards are in terms of totals; Cn is based on dissolved) 
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Remedial Scenario #2: 1996 Treatment Conditions 


1996: Gladstone treatmer 
system collected most of 
the water upstream of the 


Move upstream boundary (above Mogul) to just above South Fork 


Specify upstream boundary chemistry based on 1996 data for 
Cement Creek, just downstream of the treatment system 
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Model “Validation” 

Comparison of Treatment Simulation w/1996 Data 
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Normalized Results, Revisited 

Best Estimate Simulations 


Estimated Concentrations 

Best Estimate Simulations 


Bulkhead the Red & Bonita 
Treatment System 


A1 (total) Cd (dissolved) Fe (total) Zn (dissolved) 
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Summary 


• Neither scenario results in the attainment of chronic 
water quality standards at A72 

• Cement Creek is an appropriate focus for Fe & Zn 

• Cd: need to consider Upper Animas 

• Al: need to consider Mineral Creek 



